High energy photon-neutrino elastic scattering 
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The one-loop helicity amplitudes for the elastic scattering process 71^ 7// in the Standard 
Model are computed at high center of mass energies. A general decomposition of the amplitudes 
is utilized to investigate the validity of some of the key features of our results. In the center of 
mass, where ^/s = 2cj, the cross section grows roughly as uj^ to near the threshold for W^-boson 
production, -^s = mw- Although suppressed at low energies, we find that the elastic cross 
section exceeds the cross section for jf — > 771^ when ^/s > 13 GeV. We demonstrate that 
the scattered photons are circularly polarized and the net value of the polarization is non-zero. 
Astrophysical implications of high energy photon-neutrino scattering are discussed. 

PACS numbers: 13.15.-Hg, 14.60.Lm, 14.70.Bh, 95.30.Cq 



I. INTRODUCTION 

The scattering process "fiy has been studied 

in past using the four- Fermi theory J l| , vector boson 
theories ^, |[, the Standard Model [|, |§], and model- 
independent parameterizations At low energies, u <ti 
rue, where w is the energy of the photon in the center 
of mass and rrie is the mass of the electron, the cross 
section for the elastic scattering ji' — > is exceedingly 
small . As it is shown in the Ref . , the cross section 
for the inelastic process —^ 771^ is much larger than 
that of the elastic scattering — > jv. However, with 
increasing ui, the elastic cross section eventually exceeds 
the inelastic cross section. Here, we show that the cross 
over occurs at w ~ 7 GeV. 

Our explicit calculations, performed for high energies 
with massless electron neutrinos, show that the scale for 
71/ — > 71/ scattering is set by mw- In fact, the cross sec- 
tion is effectively independent of the electron mass, rrie 
(except near forward and backward directions). Further- 
more, since the weak interaction violates parity, we find 
a large circular polarization of the scattered photons. 

In the next section, we utilize a decomposition of 
the elastic amplitude A{s, t, u) to obtain the general 
properties of and restrictions on the helicity amplitudes 
-4aiA2 (s, u). Section III contains the numerical results 
for the complete one-loop helicity dependent differential 
and total cross sections. This is followed by a discussion 
of the production of circularly polarized photons in high 
energy 71/ collisions, and an estimate of the tem- 

perature at which, during the thermal evolution of the 
universe, photons and neutrinos decoupled. 



II. THE -yu 71/ ELASTIC AMPLITUDE 

The Lorentz-invariant amplitude for the process 'yv — > 
71/ can be expressed as |^ 

A{s,t,u) = u{p2hf,{^ + j5)u{pi)M^,al3 {ei)a{e2)0 , (1) 

where £1 and £3 ^-re the polarization vectors for the 
incoming and outgoing photons, respectively, and the 
Mandelstam variables s, t, and u are defined by s = 
-{pi + fcl)^ t = -{pi - P2f, and u = -{pi - k2f. 
Here, the 4-momenta of the incoming neutrino and pho- 
ton are pi and ki , respectively, with p2 and k2 denoting 
the corresponding outgoing momenta. The tensor Ai^ap 
can be expressed in terms of four linearly independent, 
gauge invariant tensors T^^^, i = 1, • • • , 4 |§, ^ 

+M^{s, t, u)tI^% + M^is, t, u)T^^% . (2) 

To facilitate the inclusion of the Bose symmetry, which 
requires the invariance of the amplitude under the ex- 
change of the incoming and outgoing photons, we use four 

(i) . 

tensors T^^^ which have a definite symmetry under this 
exchange. For our choice, T^^^fS symmetric, while T^J^, 



T. 



(3) 



^^p, and Tj^ij are antisymmetric. As a consequence, 
Bose symmetry requires the following relations between 
the four scalar functions Aii{s, t,u), i = 1, • • • , 4 : 

Mi{s,t,u) ^ Mi{u,t,s) , (3) 
Mj{s,t,u)= -Mj{u,t,s), j = 2,3,4:. (4) 

In the center of mass with massless neutrinos, the con- 

(i) 

tractions of the tensors T'^ct^ with the neutrino factor 
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^fj_ = u(p2)7/i(l +75)u(pi) and polarization vectors (£i)a 
and (e2)/3; result in the following helicity basis [|| 



T^MaU'^(^i)"(^2)/3 = -scos(0/2)[i(Ai+A2 + 2AiA2) 



+4SA1A2] , 

sicos(6'/2)(Ai - A2) , 
sicos(6'/2)(l - A1A2) ; 



T';jf,U^iU^*2)p = -8^cos(0/2)AiA2 



(5) 
(6) 
(7) 

(8) 



where 9 is the angle between the incoming neutrino, 
which is moving in the +z direction, and the outgoing 
neutrino, Ai = ±1 is the helicity of the incoming photon, 
and A2 = ±1 is the helicity of the outgoing photon. Of 
these contractions, Eq. (|^) is antisymmetric under the ex- 
change of Ai and A2. The other three, Eqs. (^), and 
(^, are symmetric under this exchange. The imposition 
of time reversal invariance, which implies the symmetry 
of the amplitude under the exchange of Ai and A2, means 
that the T- violating part of the helicity basis, Eq. (|^), 
must be excluded. 

Using Eqs. (0), dl), (i, (I), and (|), the hehcity am- 
plitudes A\^\2{s,t,u) can be written as 

A++{s,t,u) = 2scos{9/2)[2uMi{s,t,u) 

Sll 

_4— 7W4(s,t,u)], (9) 
A--{s,t,u) = 2scos(6'/2)[-2sA^i(s,t,u) 

Sll 

_4— 7W4(s,t,u)], (10) 

A+-{s,t,u) ^ 2scos(e'/2)[(s - u)Mi[s,t,u) 

su 

+tM3{s, t, u) + A—Miis, t, u)] , (11) 



A ^{s,t,u) = A^ {s,t,u) 



(12) 



where t — —\s{\ — z), u — — ^s(l + z), z — cos9, and 
we have assumed time reversal symmetry and, therefore, 
omitted A42{s,t,u). Notice that Eq. ( p^ is the result of 
T-invariance, not Bose symmetry. 

To include the requirements of Bose symmetry, Eqs. (|^) 
and (||), and also to exhibit the conservation of angular 
momentum in the expressions for the helicity amplitudes, 
Eqs. (|9|)-([rT|), we define the following two independent 
scalar functions 



T{s,t,u) = 4A^i(s,t,u) - 8-7^4(5, i,u) , (13) 

u 

Q{s,t,u) = Mz{s,t,u) - 2-Mi{s,t,u) 

SIX 

-h4— 7W4(s,t,u). (14) 

The interchange of s and u in the Eqs. ( [l^ ) and to- 
gether with the Bose symmetry requirements of Eqs. (^) 



and (Q), result in the following relations 

u 

J'{u,t,s) = AMi{s,t,u) + 8-i-M4is,t,u) , (15) 



Q{u,t,s) = -A^3(s,t, u) - 2-Mi{s,t,u) 



SU 

-4-^7W4(s,t,u) 



(16) 



Using Eqs. (|T|)-(|l|) in Eqs. (|)-(|ri|) results in the fol- 
lowing general expressions for the helicity amplitudes 

A++{s,t,u) = sucos{9 /2)!F{s,t,u) , (17) 
A--{s,t,u) = -s^ cos{9/2)T{u,t,s) , (18) 
A+-{s,t,u) = stcos{9/2)[g{s,t,u) - g{u,t,s)]. {19) 

In order to ensure the conservation of angular momen- 
tum in Eqs. ([l7|)-(|l9|), it is necessary to require that the 
functions T{s,t,u), T{u,t,s), and [g{s,t,u) — Q{u,t, s)] 
be non-singular in the limit u —^ (backward scattering) . 
In addition, the function [g{s, t, u) — g{u, t, s)] must also 
be non-singular in the limit i — > (forward scattering). 

From the Eqs. (p7|)-(p^, it is clear that the interchange 
of s and u, results in the following relation 



A\^\^{s,t,u) = A^\^-\^{u,t,s) 



(20) 



where, under this interchange, the factor scos(6'/2) — 
s^J—u/s becomes u^J—s/u — —s^J—u/s = — s cos(6'/2). 

We can use Eq. (^^ as a check of our calculation. To 
do this, we must express the helicity amplitudes as func- 
tions of s and u such that these functions remain well 
defined after the interchange s u. A certain amount 
of care must be exercised when performing a numerical 
check of Eq. (pO|), because it is convenient to use the fact 
that s > and u < when calculating A\^\2{s,t,u). 
The interchange of s and u would seem to move the nu- 
merical calculation into the region s < and it > in 
order to make the comparison. However, no additional 
calculation is required if one remembers that each of the 
diagrams of Fig. 1 has a counterpart in which the pho- 
tons are interchanged. If a particular scalar contribution 
to one of the direct diagrams is /(s, u) the correspond- 
ing crossed diagram will contribute f{u,s). When it is 
assumed that s > and u < 0, the result is a function 
/i(s,u) for the direct diagram and a different function 
/2(w, s) for the crossed diagram. The function fi{x, y) is 
not defined when its first variable is negative and its sec- 
ond positive, while f2{x,y) is not defined when its first 
variable is positive and its second is negative. Since the 
s <-^- u interchange essentially exchanges the direct and 
crossed contributions, we can use 

Fis,u) = 9{s)9{-u)h{s,u) + 9{~s)9{u)h{s,u) , (21) 

for the direct contribution and 

Fiu, s) = 9iu)9{~s)fiiu, s) + 9{-u)9{s)f2{u, s) . (22) 

for the crossed contribution. Clearly, F{s, u) and F(u, s) 
are well defined for both s > 0, u < and s < 0, u > 0. 
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III. DIFFERENTIAL AND TOTAL CROSS 
SECTIONS 



We calculated the diagrams of Fig. [l] in a nonlinear 
gauge such that the coupling between the photon, the 
T4^-boson and the Goldstone boson vanishes Ol. 
Using algebraic manipulation software SCHOONSCHIP |12] 
and FORM [ p3| , the diagrams were decomposed in terms 
of scalar n-point functions |l^, and then checked 
numerically with the FORTRAN codes LOOP |]l6j and ff 
IllTi. In addition, using Eqs.(|l]) and (11), we demon- 
strated that our numerical calculations for the helicity 
amplitudes, satisfy Eqs. (|l^) and (|2^). To simplify the 
calculations, we assumed that s,t,u ^ rn^. Several of 
the scalar n-point functions depended upon In(me) or 
In^(me). However, it turns out that every diagram of the 
Fig. is independent of the mg. Due to the assumption 
s, t, w 3> TTjg, in general, we do not expect that our results 
be reliable near the forward and the backward directions, 
where t —^ and u ^ 0, respectively. 

To study the degree of reliability, and also as a partial 
check of our calculated helicity non-flip amplitudes near 
the forward direction, we use the optical theorem, which 
relates the imaginary part of the non-flip amplitude for 
elastic scattering in the forward direction to the total 
cross section for the process W'^e~ as 



- - ImAxxi 
s 



= 0)=ax 



(23) 



Here, A is the helicity of the photon, and ax represents 
the total cross section for a photon of given helicity A, 
after summation over the helicities of the W-boson and 
the electron. Our explicit calculation gives 



2V2G 



pa 



+ - 



£{s, m^, ml 



m 



w 



(1 



(24) 



2V2G^ J^EfIM^(l _ !!!k)(i 



m 



w 



(1 



^2 



s 



W \2 



(25) 



where 



A(a;, y, z) = + + — 2xy — 2xz — 2yz , (26) 



n, X , / x-y + z + ^\{x,y,z)\ 
^(x,?/,z)=ln . (27) 



y + z- y/X{x,y,z) 



We retain 



the electron mass in the functions 



A(s,m^,TOg), £{s,my^,,ml) and 



£(s, TO, 



e ' 



TO 



W 



) to 



ensure that the cross sections vanish at threshold, 
-v/i = mw + rjie- In the coefficients of these functions, 
we have dropped powers of ml/m^ 



ty^. The spin aver- 
aged cross section (cr+ -I- (J-)/ 2 agrees with the result 
previously obtained by Seckel . 

Both cross sections approach the same constant when 
^/s ^ mw- Since £{s,m'^,ml) ^ ln{m'^/ml), each 
cross section contains a In(TOg) term, which our numeri- 
cal calculation oi Axxis, Q) will not obtain in the forward 
direction due to the assumption s,t,u^ m^. Neverthe- 
less, we can use the optical theorem to check the high en- 
ergy behavior of our amplitudes because the \n{m^/ml) 
terms in cr_|_ and (t_ depend very differently on to^/s; 
{my^l sY versus my^/s. Indeed, for \/s > 200 GeV, and 
A = +1, our calculations for the —ItsiA++{9 = 0)/s 
and (74- arc in good agreement with Eq. (^) , while for 
A = — 1, the cross section (T_ is almost a factor of two 

larger than the quantity —IvtvA (6* — 0)/s. However, 

at higher energies, we find an excellent agreement with 
Eq. ( [2^ ) for both A = ±1 helicities. 

Before leaving the subject of the behavior of the non- 
flip forward helicity amplitudes, it is worth noting that 
the exact value of the forward helicity amplitude Axx {s) 
can be obtained using the dispersion relation 



Axx{s) = — 



ds' f axis') , axis') 



2 5' \ S' — S S' -\- S 



(28) 



where axis) is the cross section in the channel 7P — > 
W^e^. It is not difficult to show that axis) = a^xis), 
in which case we have 



Axxis) 



„2 poa 



ds' r axis') ^ a^xjs') 

^ J{m,v+n.,)^ S' \S' - S s' + S 



29) 



This expression obeys the symmetry relation Eq. (|20| ) 
specialized to the forward direction t = 0, 



Axxis) = A-x-xi-s). 



(30) 



To leading order in s, the dispersion integral can be eval- 
uated and gives 



A+^is) ^ A__is) 



V2G 



pas 



'■w 
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(31) 



in agreement with the low energy result of Details 
will be presented elsewhere. 

In Figs.§ and || we show the differential cross sections 
for elastic scattering using 



dax-,x 
dz 



32tts 



(32) 



where, Ai and A2 are the helicities of the incoming and 
outgoing photons, respectively. Figs.^ and ^ illustrate 



4 



the identity of the amphtudcs and A y, as required 

by Eq. . They also show the vanishing of the amph- 
tudes for backward scattering, 9 = n. However, they 

do not exhibit the vanishing of the A^ or A ^ in the 

forward direction, 6 = 0. As we stated earher we do 
not, in general, expect our results be reliable near the 
forward (t — > 0) and backward {u 0) directions on a 
scale ~ mg. In our calculations, we have, for instance, 
replaced factors such as t/{t — m1) with 1. This situation 
is similar to that encountered when calculating the am- 
plitudes for quark+gluon quark-|- photon [|l9|, where 
the existence of the kinematic zeros proportional to t in 
the forward direction is obscured by the assumption that 
the quarks are massless. 

The total cross sections, for helicities Ai and A2, are 
given by 



1 



327rs 



dz\A\^> 



(33) 



and are plotted in Fig.Q. Shown in dots is the cross 
section for helicity flip scattering, which can be seen to 
be much smaller than the cross sections for the helicity 
non-flip scattering. This feature seems not to be a con- 
sequence of any symmetry, but it is reminiscent of the 
low energy case Ufs <C me), where the helicity flip am- 
plitudes vanish Q. In Fig.||, we show the total cross 
section for an unpolarized initial photon. 



_ 1 . 



+ (T+_ +cr-+) 



(34) 



This figure illustrates the roughly behavior of the total 
cross section to near the threshold for the VF-boson pro- 
duction. A fit to the points in Fig. ^, for <C w <C m^i/, 
yields 



yi. = 6.7 X 10 



-33 



pb. 



(35) 



where lo = is the energy of a photon (or a neutrino). 

In the Ref . |2y] , it is shown that the cross section for the 
process — > jjiy, in the range of energies rrie <ti uj <ti 
mw, can be written as 



1.74 X 10- 



-16 



pb. 



(36) 



Comparison of the Eq. ( |35| ) with the Eq. (^) , shows that 
the two cross sections are equal for u = 1.27 x lO'^rrie 
or ^/s = 13 GeV. Therefore, at sufficiently high energies, 
the process 71^ dominates the process "fv — s- jjiy. 



IV. DISCUSSION AND CONCLUSIONS 



when ^/s ^13 GeV. Since the weak interaction violates 
parity, the final photons in both processes acquire circu- 
lar polarization. In the case of inelastic scattering, it is 
shown in the Refs. ||^ and ||2^ that the circular polariza- 
tion is of order 20 — 30% for center of mass energies less 
than lOOme. 

To assess the degree of circular polarization of the final 
photon in the elastic scattering, we define the polariza- 
tion P as 



or, smce (t_|__ 



(7 + cr+- + CT-H 

0"-+ <C cr , 



a 



(74 



(37) 



(38) 



In Fig.^, we have plotted the polarization P as a func- 
tion of the center of mass energy ^/s. As is clear from this 
graph, for wide range of energies above 1 GeV, the polar- 
ization is of order 20 — 30%, while for energies around 200 
GeV it can reach over 60%. To find P for center of mass 
energies y/s <C 2me, we can use the helicity amplitudes 
in Eq. (4) of Ref. |], to obtain P = 1/3. 

The angular dependence of the the final photon's po- 
larization can be obtained from P{z) defined as 



P(z) 



da / dz + dCT-i / dz — da |_ / dz — da++ / dz 



da / dz + da^ /dz + da |- /dz + d(T++ / dz ' 

(39) 

or, since da^ /dz = da \-/dz, 



P(z) = 



da /dz — da++/dz 



da--/dz + 2da^ /dz + da+^ / dz 



(40) 



The polarization P{z) is plotted in Fig.^ as a function 
of z = cos 61. In this figure, the solid line is effec- 
tively unchanged for the range of center of mass ener- 
gies 1 GcV ^ ^/s < 30 GeV. The dashed curve shows 
that the forward amplitudes A++{s) and A--{s) are 
not equal above the threshold for W production. This 
is consistent with the dispersion relation Eq. (p9|). Also 
included in this figure is the polarization for the case 
^ 2me, which is based on the Eq. (5) of the Ref. ||] 

(p(z) = -v{e)). 

It is worth noting that both the behavior of the 
elastic cross section and the angular dependence of P{z) 
obtained in the low energy limit |^ persist to energies 
of order mw- This means that a low energy effective 
interaction of the form lEl| 



(41) 



where T^^ 



and T;^^ 



As shown above, the cross section for elastic scattering 
is larger than the cross section for the inelastic scattering energies. 



are the symmetrical energy- 
momentum tensors of the neutrino and the photon, gives 
an accurate description of elastic scattering to quite high 
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Finally, to investigate the importance of the reaction 
71/ —^ in cosmology, we define a^, by 



f 



n^ct 



(42) 



where n^, is the neutrino number density (number of neu- 
trinos per unit volume), c is the speed of light, and t is 
the expansion time of the universe. Taking Ui, = 56 cm~^ 
and t = 15 X 10^ years, we find the present value of ct^ 
to be 



= 1.26 X 10^ pb. 



(43) 



The mean number of collisions between a photon and 
relic neutrinos |2^ is of order (J^u^-tv / (^v From the Fig.^ 
and Eq. ( ^ ) we find that, regardless of the center of mass 
energy, a^u^^^/o'v ^ 10~^. Therefore, the process 
71/ effectively ceased to occur early in the evolution of 
the universe. 

To estimate the time (or the temperature) at which 
the decoupling of photons and neutrinos in the process 
71/ jh' took place, we must determine when the value 
of the ratio a-^^-^-^^l Ui, is of order 1, 



tvl^v ^ 1 



(44) 



To do this, consider the product a{p^,p^) v^^, where 
a{pj,p^) is the cross section for the scattering of a pho- 
ton of momentum p-y with a neutrino of momentum p^, 
and is the magnitude of their relative velocity. We 
define the average value of this product by 



{a{pj,p^)v^^) 



J drij dnva{pj,p^)vji, 



where 



dn. 



1 



2d^p 



7 



drii^ 



^ (27r)3 e^^/T _ 
1 d^p^ 



(27r)3 e-E./T + 1 ■ 



(45) 

(46) 
(47) 



Here, and E^, are the energies of the photon and neu- 
trino, respectively, and the factor 2 in the Eq. (^) is 
due to the number of helicity states of the photon. Af- 
ter performing the integration in the denominator of the 
Eq. (|45|), we obtain 



{a{p^,p^)v^^) 



1 



dn~fdn^a(j}ry,pi,)v~f^ , (48) 



where and n^, the number densities for the photons 
and neutrinos, respectively, are 



1 



(27r)3 
1 

(2^ 



2d^p^ 



eE-,/T _ 1 



2C(3)T3 

3C(3)r3 

47r2 



(49) 
(50) 



Using the invariance of a{py,p^) v-fy E^ E^, the prod- 
uct a{p.y,pi,) can be expressed in terms of the cross 
section (Tcm in the center of mass as 



0'(P7,Pi') = CTCM 



P2 

E^Ei, 



(51) 



where ctcm — u^i^^^i, is given by Eq. ( |35| ) with uo = Eqm ■ 
The center of mass energy Ecu for a photon (or a neu- 
trino) in terms of E^, E^, and 0^^, the angle between p^ 
and Pi, , is 



(52) 



EcM = \/E^E„ sm{e^y/2) . 
Therefore, Eqs. (||), (|l|), and (g) give 

2E^E^ 

a{pj,py)vjy = 6.7 X 10"^^ — sm^{0^y/2) ph. 

1-1 . . ^^^^ 

Using this result in Eq. ( [48| ) and performing the integra- 
tion, we find 



{a{pj,py)vji,) = 6.7 X 10 



-33 



124 



^12 



59535 [C(3)]2 



^6 

— Pb. (54) 



The average number of collisions that a photon makes 
with neutrinos through the process during the 

time t is 



(55) 



Similarly, the average number of collisions that a neutrino 
makes with photons through the process "fv — > during 
the time t is 



Afy = {a{p^,py)v^i,) n^t . 



(56) 



For the duration time, we take an expansion time, which 
is given by [EJ 



2 X 102" 



(57) 



where T is the temperature in Kelvin. To estimate the 
temperature at which the thermal decoupling of the pho- 
tons and neutrinos took place, we use the criterion 



max[A/'T,,7V^] ^ 1 . 



(58) 



It is clear from Eqs. (^) and (|50| ) that > riy, and from 
Eqs. (11) and (|6|) that < Mu- Therefore, Eq. (H) 
gives 



Using Eqs. (p 



AC - 1. 
5|), and 



3), we find 



AC = 2.5 X lO^^^T'' 



(59) 



(60) 



Thus, from Eq. ( |59| ) we find the decoupling temperature 
T - 1.2 X 10" K, or T ~ 1 GeV. This translates into 
an expansion time of about 1.4 x 10^^ s. The tempera- 
ture T ^ \ GeV corresponds to a center of mass energy 
^/s ^ 2 GeV, which is well within the range of validity 
of Eq. (PI). 
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[22] See Chapter 6 of P. J. E. Peebles, Principles of Phys- 
ical Cosmology (Princeton University Press, Prince- 
ton, NJ, 1993). 

[23] Strickly speaking, the coefficient 2 in this equation 
should be replaced by a factor less than 1, since the 
temperature at which we are interested is of the or- 
der 10^'^ K. For a discussion, see Chapter 15 of S. 
Weinberg, Gravitation and Cosmology (John Wiley 
& Sons, NY, 1972). 
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FIG. 2: The helicity dependent differential cross sections for 71/ — > 72/ are shown for ^/s = 20 GeV. The solid line is 
da /dz, the dashed line is daj^j^/dz, and the dotted line is daj^ /dz = da-+/dz. 
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FIG. 3: Same as Fig. | with ^ 200 GeV. 




0.1 1.0 10.0 10 100 1000 10000 



Vs" (GeV) Vs" (GeV) 

FIG. 4: The hehcity dependent total cross sections for — > 71/ are shown. The sohd hne is (T__, the dashed hne is 
and the dotted hne is (T-| = cr ^. 
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FIG. 5: The total cross section a-^y^-fi, for unpolarized photons is shown. 
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FIG. 7: The circular polarization P(z) of the final photon in the process "fv — > 71/, as defined in the Eq. (|4C|), is shown. 
The solid and the dashed lines are polarization for the center of mass energies of 20 GeV and 200 GeV, respectively, 
while the dotted line is for ^/s ^ 2me, which is taken from the Ref. 



